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The problem of continuum percolation in dispersions of rods is reformulated in terms of weighted 
random geometric graphs. Nodes (or sites or vertices) in the graph represent spatial locations oc¬ 
cupied by the centers of the rods. The probability that an edge (or link) connects any randomly 
selected pair of nodes depends upon the rod volume fraction as well as the distribution over their 
sizes and shapes, and also upon quantities that characterize their state of dispersion (such as the 
orientational distribution function). We employ the observation that contributions from closed loops 
of connected rods are negligible in the limit of large aspect ratios to obtain percolation thresholds 
that are fully equivalent to those calculated within the second-virial approximation of the con¬ 
nectedness Ornstein-Zernike equation. Our formulation can account for effects due to interactions 
between the rods, and many-body features can be partially addressed by suitable choices for the 
edge probabilities. 


I. INTRODUCTION 

Dispersions of rod-like particles in the continuum are 
prototypical representations of numerous and varied sys¬ 
tems of diverse and widespread technological interest, 
including liquid crystals, fiber-reinforced materials, and 
polymers filled with conducting nanofibers and nan¬ 
otubes. Of particular interest are percolation phenomena 
associated with long-range connectedness of rod-like 
particles, as these can dramatically affect the strength of 
Hber-reinforced materials and the electrical conductivity 
of nanocomposites. 

Traditionally the volume fraction at which an infinite 
cluster of connected particles first arises (referred to as 
the “percolation threshold”) has been studied theoret¬ 
ically by means of integral equation approaches based 
upon the connectedness Ornstein-Zernike equation M- 
For the case of elongated rods that are distributed uni¬ 
formly and oriented randomly in the matrix, the predic¬ 
tion that the percolation threshold depends in an approx¬ 
imately inverse manner upon the aspect ratio, has been 
confirmed by computer simulations for both penetrable 
and hard-core-soft-shell rods 

In recent years, the impact of polydispersity in the 
particle size (the fact that in real-life situations, parti¬ 
cles usually display a distribution over a range of sizes 
and are not all of uniform dimensions) upon the perco¬ 
lation threshold has been examined using both (i) inte¬ 
gral equation methods ii , and (ii) a mapping between 
continuum percolation and percolation on an appropri¬ 
ately modified, tree-like, Bethe lattice [l^. The latter 
formalism applied excluded volume arguments [l^ to es¬ 
timate the average number of inter-particle contacts as a 
function of volume fraction, and related these to nearest- 
neighbor occupancy rates on a Bethe lattice with suit¬ 
ably chosen vertex degrees (co-ordination numbers). For 
cases in which the particles are assumed to have polydis¬ 
persity in their lengths alone and uniform diameters and 


widths, both of these approaches show that the leading- 
order term determining the percolation threshold varies 
inversely with the weight-averaged aspect ratio of the 
particles, as confirmed recently by Monte Carlo simula¬ 
tions [I^HjliSl- The original implementations of the lat¬ 
tice analogy neglected fluctuations in the expected num¬ 
ber of contacts between particles, and led to findings that 
differed from those obtained from the integral equation 
methodology for non-leading-order terms in the expres¬ 
sion for the percolation threshold. 

In this work, the continuum percolation problem for 
rods is reformulated in terms of random geometric graphs 
with appropriately specified degree distributions, a vehi¬ 
cle that is usually employed to describe connectedness in 
assemblies of isotropic, spherical objects (Tsl - f^ . In the 
present formulation, nodes (or sites or vertices) represent 
rods that are distributed randomly throughout the ma¬ 
trix, and the probability that a randomly chosen pair of 
nodes is connected by an edge is weighted by the rela¬ 
tive orientation of the rods. It is shown that the node¬ 
degree distribution automatically includes fluctuations in 
the number of inter-particle contacts, and that closed 
loops of intersecting rods can be neglected in the limit 
of large aspect ratios. Results for the percolation thresh¬ 
old obtained from this formulation are in full agreement 
with those previously obtained using the meth¬ 

ods of integral equation theory. In addition, it is shown 
that many-body effects for interacting rods and/or non¬ 
equilibrium features in the distribution of particles can 
be partially accounted for by suitable choices for the edge 
probabilities. 

We begin this account by first examining the case of 
fully penetrable, monodisperse rods in Sec. m and use 
this relatively simple situation to develop our notation. 
Percolation in systems of polydisperse but interpenetra- 
ble rods, and percolation in systems of rods with interac¬ 
tions (such as pairwise excluded volume), are discussed 
in Secs. m and HVl respectively, and Sec. |V] presents a 
summary of and conclusions from this study. 
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II. MONODISPERSE PENETRABLE RODS 

A. Mapping onto a random graph problem: 

Degree distributions 

We start by considering a set of N rod-like parti¬ 
cles, the centers of which are randomly located within 
a three-dimensional region of volume V. In the interests 
of achieving a clear presentation, and in order to establish 
our notation, we start by considering hrst and separately 
the case of monodisperse rods, which we model as iden¬ 
tical penetrable cylinders of length L and diameter A. 

The orientation of each rod is specified by a unit vec¬ 
tor, ft, that is aligned parallel to the axis of the cylinder. 
We assume that each unit vector u may take on any one 
of M„ < N possible orientation Ui (with i = 1 , 2 ,..., Mu 
) that are drawn from a prescribed distribution function. 
A pair of cylinders with orientations Ui and Uj and cen¬ 
tered at the locations r and r' are said to be connected to 
each other if the (penetrable) boundaries of the particles 
overlap each other. We next construct a network with N 
points (or vertices or nodes) by associating a node with 
the center of each particle. Edges (or links) connecting 
pairs of nodes are assigned based upon whether or not a 
given pair of particles (or nodes) are connected to each 
other based upon the foregoing criterion for inter-particle 
overlap. This A-point network may be regarded as a 
weighted random geometric graph, in which the proba¬ 
bility of edge formation between any pair of randomly 
selected points (nodes) depends upon the relative posi¬ 
tions and orientations of the cylinders that are associated 
with the selected nodes. We next introduce the connect¬ 
edness function = f{R,Ui,Uj), where R = f — r' is 
the displacement vector between the centers of the pair 
of cylinders, defined such that fij is equal to unity if the 
pair of cylinders overlap and is otherwise equal to zero. 
The probability that the nodes representing particles of 
orientations Ui and Uj are directly linked by an edge is 
then given by the integral of fij over R\ 

Vij = — J' dRf^R^Ui^Uj'). ( 1 ) 

The displacement vector can be expressed as i? = zui + 
z'uj+X{uiXUj)/\uiXUj\ [2^, where A is the component of 
R along the direction of the shortest line connecting i and 
j, and 2 and z' are components along the directions of Ui 
and Uj. In the limit L 3> A, the connectedness function 
is expressed in terms of z, z', and A as simply fij = 
e{L /2 - \z\) 9 iL /2 - \z'\)e{A - |A|) glU , where 6 is the 
Heaviside step function. Using dR = | siny^-Idzdz'dA, 
where 7^ = 'y{ui,Uj) is the angle between the directions 
of Ui and Uj, Eq. m reduces to Vij = V^^fV, where 

U,f^ = 2L2A|sin7„-| (2) 

is the leading-order contribution to the excluded volume 
[2^. The dependence of upon the relative orienta¬ 
tion of the pair of cylinders can be thought of as providing 


a weight factor for the edge probability that is maximum 
when the cylinders are orthogonal and vanishes for par¬ 
allel configurations. Note that the vanishing of Vij when 
= 0 and tt is an artifact of the L ^ A limit, which 
has negligible impact on the connectivity of slender rods 
with isotropic orientational distributions. 

In order to proceed further, we must introduce the 
degree distribution function that characterizes the con¬ 
nectedness properties of the network. We construct the 
probability pi{ki^ui\ k2, M2; • ■ •) that a randomly selected 
cylinder i with orientation Ui is connected to fci cylin¬ 
ders with orientation mi, k2 cylinders with orientation 
U2, and so on. Since we have assumed that the orienta¬ 
tions of different cylinders are independently distributed, 
Pi{ki, ui; ^2, M2; • ■ •) is a product of binomial distributions 
Pij {kj), each of which gives the probability that the cylin¬ 
der of type i is connected to precisely kj cylinders of type 

j- 

Pj(A:i,mi; fc2,M2; • ■ •) = (3) 

j 

where 

( 4 ) 

where Nj is the number of cylinders with orientation vec¬ 
tor Uj, and dij is the Kronecker symbol. We next take 
the limit U —>■ 00 in such a manner that the number den¬ 
sity Pi = Ni/V remains finite for all i. Introducing the 
number fraction Xi = Ni/N and the total number den¬ 
sity p = N/V, Eq. ([ 4 ]) reduces to a Poisson distribution 
in this limit: 

kj 

P^jikj) = ( 5 ) 

Kj. 

where 

Zij = XjpVfj^^ = 2xjpL^A \sinyy I, (6) 

is the average number of cylinders of type j that overlap 
a given cylinder of type i. The distribution function in 
Eq. ([ 5 ]) accounts for fluctuations in the number of con¬ 
tacts that may exist between pairs of cylinders of types 
i and j, which is an issue that was neglected in a prior 
effort to model percolation by rods within a lattice-based 
approach [T5j| . 

Knowledge of the degree distribution of Eq. is suffi¬ 
cient to fully ascertain the connectivity properties of the 
system, provided that we consider the limit of slender 
rods (that is, L ^ A) and assume that the orientations of 
the rods are distributed isotropically. This is so because 
under these conditions the probability of Ending closed 
loops in finite clusters of isotropically oriented connected 
cylinders scales as A/L (as shown in the Appendix) and 
thus vanishes in the limit L/A —>■ 00. Consequently, the 
random geometric graph constructed as described above 
has a tree-like, dendritic, structure, where the number of 
branches attached to a randomly selected node follows 
the distribution function given by Eqs. m and ®. 
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B. Calculation of the percolation threshold as 


In this section, we calculate the percolation thresh¬ 
old by following a method that takes advantage of the 
tree-like structure of the graph developed in Sec. Ill Al 
We consider a randomly selected edge that connects a 
node of type i with a node of type j, and we follow that 
edge from i to j. The node of type j that we arrive 
at by following that edge will be ki times more likely 
to have degree ki than degree one in terms of its direct 
links to nodes of type i. The degree distribution for the 
node of type j that is arrived at in this way, denoted 
qji{ki^ui-,k 2 ,U 2 ] ■ ■ ■), must therefore be proportional to 
kiPj{ki,ui;k 2 ,U 2 : ■ ■ ■) [25l - [^ . Using Eqs. ((3|) and ([5]) 
and requiring that qji{ki,ui \ k 2 ,U 2 ] ■ ■ ■) be correctly nor¬ 
malized leads to 

qp{ki,ui]k2,U2\. ■ ■) = WPjiih) 

= Pji{ki-l)\Ypji{ki). (7) 

Given the dendritic structure of the network, we can 
write the probability Qji that none of the remaining other 
edges leading out of the node of type j leads to the the 
giant component as follows: 

9A(fccWi;fc2,M2;...) (8) 

ki,k2,... Tfi 

where the Kronecker symbol Smi accounts for the fact 
that there are ki — 1 edges with nodes of type i other 
than the initially selected one that was followed to arrive 
at the node of type j. Equations 0 and dH) lead to 


1)]^ 


L k 


^Pji{k)Q 


= exp 


ZjiiQij 1 ) 


(9) 


from which we see that Qji depends only on the first in¬ 
dex: Qji = Qj. It should be noted that Eq. always 
admits the trivial solution Qj = 1, which corresponds to 
the absence of infinite, connected clusters. The existence 
of non-trivial solutions such that Qj < 1 indicates the 
presence of an infinite cluster (giant component) com¬ 
prising a finite fraction of nodes, and that spans the en¬ 
tire system [ij. The percolation threshold is identified 
by evaluating the smallest fraction of nodes for which 
such a non-trivial solution of Eq. 0 exists. We thus set 
Qj = 1 — Ej, with Ej <IC 1, and expand Eq. 0 to first 
order in Ej , which yields 


Sj =^Zji€i = 2pL‘^A^xi\siiijji\ei, (10) 

i i 

where we have used Eq. 0. Passing to a continuum 
representation of the orientational distribution by intro¬ 
ducing p(u) = '^iXi5(u — Ui) , Eq. (flUl) can be written 


E{u)=2pL^A J d'u'p('u')| sin[ 7 ('u,'u')]|e:('u'). (11) 

Eor the case of an isotropic orientational distribution of 

the rods, for which p(u) = l/47r, performing the integral 
over u on both sides of Eq. (HD, we find the following 
result for the percolation threshold: 

ttP^L^AI2 = 1 , ( 12 ) 

where we have used (47r)“^ J du| sin[ 7 (u, u')]| = 7r/4 
and Pc denotes the number density at the percolation 
threshold. Introducing the dimensionless critical density 
Tic = Pc^L/4:, Eq. (IT^ can be rewritten as 

= (13) 

which coincides with the percolation threshold of slen¬ 
der rods obtained from connectedness percolation theory 
within the second-virial approximation [^. 

We have obtained Eq. dT51) by finding the density of 
rods such that Eq. (O has a non-trivial solution (other 
than the trivial solution identically equal to unity). Ow¬ 
ing to the absence of correlations between penetrable 
rods, we could, instead, have followed a different proce¬ 
dure by initially choosing a density p' for the rods that is 
(by construction) large enough that an infinite connected 
cluster is certain to exist, and subsequently randomly 
deleting a given fraction q of rods until the damaged net¬ 
work just barely ceases to remain infinitely connected 
(the giant component disappears) [l^. This alternative 
procedure is equivalent to rewriting Eq. ([8]) as 

Q]i = q+{^-q) X] 9a(^i>'“i;^2,u2;...)P(5^"‘^'”% 

(14) 

where qji{ki,ui] ^ 2 , U 2 \ ■ ■ ■) is the distribution of Eq. ([7]) 
for rod density fixed at p'. By repeating the same steps 
of the above analysis, we find that the critical fraction qc 
of deleted nodes satisfies 

(1 - qc)Trp'L^A/2 = 1, (15) 

so that by defining 77 c = (1 — qc)p' tvA^L/A we obtain a 
result identical to that in Eq. m- We note that the 
above procedure of randomly deleting nodes from a pre¬ 
existing gian t component is equivalent to the one followed 
in Refs. [l5j to establish a mapping between continuum 
percolation of rods and percolation on a modified Bethe 
lattice. In the present approach, however, fluctuations in 
the number of contacts are fully taken into account, as is 
apparent from Eq. 0. 


III. POLYDISPERSE PENETRABLE RODS 

Having established our notation and formalism for the 
case of monodisperse rods, we next examine situations 
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in which the cylinders have a distribution of lengths and 
diameters. We start by adopting a discrete representa¬ 
tion, such that the length and diameter of each cylin¬ 
der may assume the values Li^ and respectively, 
with = 1, 2,... and ja = Ij 2,..., with probabili¬ 
ties denoted The three-component index vector 

i = (li, lA, iu) is introduced to identify the length, diam¬ 
eter, and orientation of each cylinder. With this notation, 
a straightforward generalization of Eq. o shows that the 
probability that a pair of such randomly located and pen¬ 
etrable cylinders are connected is given by uy = 
where 


For an isotropic orientational distribution for the rods, 
averaging over the angles leads to 

II ,1a 

+ X! (22) 

II ,1a 

where Introducing auxiliary vari¬ 
ables defined as: Si = Ja ^2 = 

^JL,JA ^JLjALjL^jA£jL,jA, Eq. ^ cau be rewritten as 


I sin7,„,,„ I (16) 

is the excluded volume for a pair of cylinders 
with lengths, diameters, and orientations given by 
Li ^, Ai^, Ui^ and Lj ^, Aj^, Uj^, respectively. (The 
derivation of Eq. (1161) implicitly assumes that the lengths 
of the cylinders are much larger than their diameters.) 
An analysis very similar to that presented in Sec. Ill Al 
shows that the probability Pi(fcj, j; k\, 1; ...) that a node 
of type i is connected to precisely fcj nodes of type j, k\ 
nodes of type 1 , and so forth, is 


£1 = Jp(T^A)£i + \p{L^)e2, 

82 = + Jp(T^A)£2, (23) 


where (L^A™) = a;ji,jA(^iL)^(AA)™ for "i = 

0, 1, 2. The condition that the linear and homogeneous 
system of Eq. (12311 has a non-trivial solution may be ex¬ 
pressed in terms of the vanishing of the appropriate de¬ 
terminant formed from the coefRcients, that is. 


(Jp,(L 2A) - l)' - (Jp,)' {L^){L^A^) = 0, (24) 


Pi(fcj,j;fci,l;...) =npij(fcj), 
j 


(17) 


from which we find that the critical threshold pc = 
■KpciLA^)/A is given by 


with 


kj 

(18) 

where 

= “^XjLjA^juPLiLLjL I sin7,„7„| (19) 

is the average number of edges that an i node forms with 
nodes of type j. Within the ansatz of tree-like connected¬ 
ness for the rod network, analysis that follows the steps 
of Sec. Ill Bl leads to the followine analoe to Eo. that 
now includes the index vectors: 


(EA^) 

(L2A) + v/(T2)(T2A2)’ 


(25) 


which coincides with the result obtained from contin¬ 
uum percolation theory based upon the connectedness 
Ornstein-Zernike equation Q. Note that for cases in 
which the rods are assumed to have polydispersity in 
their lengths alone and identical diameters, Eq. (1251) re¬ 
duces to Pc = where Lu, = (L^)/(L) is the weight 

average of the rod lengths. It is straightforward to show 
that Eq. (1^ can also be obtained by randomly deleting 
a fraction qc of nodes from an infinitely connected cluster 
until the giant component ceases to exist, as was done in 
Sec. Ill Bl for the case of monodisperse rods. 


IV. INTERACTING RODS 


Qji = exp 


^^2ji(Qij 1) 


( 20 ) 


Setting Qji = Qj = 1 — ej, we find the following for 
£j <C I: 


1 

(21) 


In actual real-life systems containing rod-like particles, 
such as nanofibers or nanotubes dispersed in polymeric 
matrices, hard-core potentials that are reflected in steric 
and excluded-volume effects prohibit the intersection of 
pairs of particles. This consideration limits the practi¬ 
cal utility of the model of ideal, fully penetrable rods. 
Furthermore, surface functionalization or addition of de- 
pletants can induce effective attractive or repulsive forces 
between the rods that strongly influence the onset of per¬ 
colation [n|, [ 29 I - I 31 II . In order to illustrate the effects of 
rod interactions upon connectedness properties within 
our present formalism of random geometric graphs, we 
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consider for simplicity impenetrable monodisperse cylin¬ 
ders with the hard-core diameter D, surrounded by a 
cylindrical penetrable shells of thickness (A — D)/2 > 0. 
As before, we work within the limit where the rods can be 
treated as being slender or of high aspect ratio, namely, 
that L A, D. For homogeneous dispersions of rods 
the probability that a pair of cylinders of orientations Ui 
and Uj are directly linked by an edge is given by 

Vij = ^ J dRf{R,u^,Uj)g2{R,Ui,Uj), (26) 

where g 2 {R,Ui,Uj) is the pair distribution function de¬ 
fined such that V~^g 2 {R,Ui,Uj)dR is the probability of 
finding a rod of orientation Uj within the volume ele¬ 
ment dR centered about the position of a randomly se¬ 
lected rod of orientation Ui [^. For ideal, penetrable 
cylinders (that is, for D = 0), g 2 {R,Ui,Uj) = 1, from 
which we recover Eq. O- For impenetrable cylinders, 
g 2 {R,Ui,Uj) = 0 whenever the hard cores of the two 
cylinders overlap. In the slender rod limit and for small 
rod densities, we take the pair distribution function to 
have the following approximate form: 

g2{R,d,,u,)^e{\X\-D), (27) 


where A is the component of R along the direction of the 
shortest line connecting i and j. The form of g 2 {R, Ui, Uj) 
in Eq. (|27)l forbids interpenetration of the cores of the 
cylinders and assumes that for |A| > D the rods are 
completely uncorrelated. Using the coordinate system 
introduced in Sec. Ill A1 Eqs. (E51) and (E71) yield 


_ 27.2 (A _ £))| sinqyl 

- y 


(28) 


which corresponds to the L ^ A, D limit of the excluded 
volume in units of V for a pair of hard rods with penetra¬ 
ble shells. Following the same steps described in Sec. m 
and assuming that the node degree distribution can still 
be expressed as a product of Poisson distributions as in 
Eqs. ([31) and (|4|), we find that the critical volume frac¬ 
tion for the cores of the rods = t^PcLD'^/A at which 
the system first percolates is 




1 D2 

2 L(A - D) ’ 


(29) 


which coincides with the percolation threshold calculated 
by solving the connectedness Ornstein-Zernike equations 
for hard-core-soft-shell rods in the second-virial approx¬ 
imation Analysis performed along lines that are 

similar to those followed in this section and in Sec. uni 
reveals that, for isotropically oriented polydisperse rods 
with impenetrable hard cores, our formalism yields a re¬ 
sult that is identical to that obtained from the connect¬ 
edness Ornstein-Zernike equations in Ref. 

Our description of the impact upon the percolation 
threshold of addressing inter-particle interactions can be 


generalized by rewriting Eq. dUl) as follows: 

g 2 {R,Ui,Uj) = x(|A|,(())6»(|A| - D), 


(30) 


where x(|A|,(^) allows for correlations between hard rods 
even for |A| > D (for illustrative purposes, the following 
analysis neglects possible dependencies of x on the ori¬ 
entation vectors Ui and Uj). In principle, xd-^lj^^) may 
also describe cases in which impenetrable rods interact 
through an additional short-range attractive or repulsive 
potential and/or nonequilibrium features of the particle 
distribution. It is straightforward to show that Eqs. (IMl) 
and (1501) lead to the following implicit equation for the 
critical volume fraction: 


Xa(/'c)/'c 


1 £>2 
2L{A-Dy 


(31) 


where xa(/') = (A — D)~^ dXx{X,(j)). Although 

knowledge of (pc depends on the specific functional form 
assumed for xa (</>), Eq. (1511) may nevertheless be used 
to predict the main effects of correlations between rods 
on the percolation threshold. For example, for hard 
rods with sufficiently small values of A — £, we can ap¬ 
proximate xa{4>) by the contact value of the pair dis¬ 
tribution function. Using either the Carnahan-Starling 
approximation for hard spheres at the same volume 
fraction XA(fP) = (1 — ^/2)/(l — (/)^, as proposed in 
Refs. [^,1^, or the Lee-Parsons approximation xa{4>) = 
(1 — 3(//4)/(l — (/)2, as suggested in Ref. [13 , we see from 
Eq. (I5T|) and from the observation that xa(</>) > 1 that 
the percolation threshold is systematically reduced when 
compared to Eq. (1551) as many-body effects are taken into 
account. 


V. DISCUSSION AND CONCLUSIONS 


By extending the formalism of geometric random 
graphs to the case of random dispersions of penetra¬ 
ble rods with isotropic orientations, we have obtained 
percolation thresholds that are fully equivalent to those 
calculated within the second-virial approximation of the 
connectedness Ornstein-Zernike equation i 13. This 
result rests on the observation that contributions from 
closed loops of connected rods can be safely neglected 
for rods with large aspect ratios, so that graph compo¬ 
nents have a tree-like, dendritic structure that allows for 
a straightforward analytic solution for locating the con¬ 
ditions that correspond to the onset of percolation. Our 
analysis of the case in which rods have an impenetra¬ 
ble hard-core has evidenced that a similar equivalence 
holds between integral equation theory and random ge¬ 
ometric graphs. In particular, the second-virial approx¬ 
imation of the connectedness Ornstein-Zernike equation 
corresponds to taking node-degree distributions that are 
products of Poisson distributions for each orientation of 
the rods. 
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We conclude by pointing out that, in contrast to inte¬ 
gral equation approaches to continuum percolation ii, 
connectedness in random geometric graphs is not re¬ 
stricted to equilibrated and random homogeneous dis¬ 
tributions of particles. The present formalism is (at least 
in principle) also applicable to more general systems, in 
which kinetic and non-equilibrium effects that may re¬ 
flect the processing history of the material modulate the 
frequency with which contacts are formed between par- 
ticle pairs, and thereby the percolation threshold 


Appendix A: The n-cycle coefficient 

Here we estimate the probability of finding a closed 
loop of n intersecting rods by considering the n-cycle co¬ 
efficient, which is defined as the conditional probability 
that two rods are connected given that they are the end 
nodes of an n-chain of connected rods. For penetrable 
rods of identical length L and diameter A, the n-cycle 
coefficient is given by 


/ n 

dridUip{Ui)f 12/23 • • • /n-l,n/n,l 

_ 

/ n 

dridUip{ui)f 12 f23. ■ ■ ■ fn-l,r 


(Al) 


where fij = f{rij,Ui,Uj) is the connectedness function 
introduced in Sec. Ill Al for a pair of penetrable rods 

I 


with orientations Ui and Uj with their centers sepa¬ 
rated by fij = fi — fj. In writing Eq. (lAll) we use 
the orientational distribution function p{u), introduced 
in Sec. IIIBI which satisfies the normalization condi¬ 
tion J dup(u) = 1. Introducing the Fourier transform 


f{q,Ui,Uj) = f dff{f,Ui,Uj)e , Eq. (lAll) reduces to 

/ d/ f ^ 

( 2 ^ J fiq,Un,Ui) 


Cn — 


dUip{Ut)f{0, Ui,U2)--- /(O, Un-l,Un) 


(A2) 

In the limit of slender rods (L ^ A), we calculate the 
leading-order contribution to /(g, Ui, Uj) by using the co¬ 
ordinate system introduced in Sec. Ill Al 


f{q,Ui,Uj) =2L^A| sin7y |jo(g • UiL/2)jo{q- UjL/2) 

X jo[Aq- {ui X Uj)/\ui X uj|], (A3) 


where jo{x) = sin(a::)/x. We next introduce the dimen¬ 
sionless momentum defined by y = gT/2 and approxi¬ 
mate the last factor in Eq. (IA3I) as jo[{2A/L)y ■ {ui x 
Uj)/\ui X Uj\] « 1 for {2A/L)y\aij\ < 1 and zero other¬ 
wise, where = y ■ (ui x Uj)/\ui x Uj\. In this way, 
Eq. (IA2I) can be rewritten as 

Cn ~ (Ad) 


where 


J J Wduip{ui)9(^ - ^amax2/^ io{y-ui?' • • •io(27- sin7i2| • • • |sin7„i| 

P n ■ 

/ J^duip(uj)| sin7i2| • • • | sin7„_7„| 


(A5) 


where amax is the maximum value that can be assumed 
by any of the terms |ai 2 |, |q; 23 |, ■ • ■, \ani\- For the case 
of an isotropic orientational distribution of the rods, the 
denominator of Eq. (IA5p reduces to (7r/4)"“^. Noting 
that I sin 7 y I < 1, we find the following upper bound for 
Cn- 


Cn < 


— jj-n+2 

22n+l 

TTlt+l 


/ 


_ , 2A 
dy9{ 1 - —Ctma^y 




rL/ 2 A 


dy 


dtjoitf 


uo 


(A6) 


where we have set amax ~ 1 in the last integral. In the 
limit of large aspect ratios, the above integral becomes 


proportional to ln(L/2A) for n = 3, while it is indepen¬ 
dent of L/A for n > 3. We thus obtain that the leading 
order contribution to Eq. (IA4I) scales as 


Cn oc 



for n = 3, 
for n > 3, 


(A7) 


from which we recover the estimate for n = 3 given in 
Ref. [i^. Equation (IA7I) shows that, for rods with suffi¬ 
ciently large aspect ratios, the impact of accounting for 
such closed loops in calculating the percolation thresh¬ 
old is likely to be minor, and to diminish with increasing 
aspect ratios. 




















7 


[1] D. Stauffer and A. Aharony, Introduction to Percolation 
Theory (Taylor and Francis, London, 1994). 

[2] S. Torquato, Random Heterogeneous Materials: Mi¬ 
crostructure and Macroscopic Properties (Springer, New 
York, 2002). 

[3] A. L. R. Bug, S. A. Safran, and I. Webman, Phys. Rev. 
Lett. 54, 1412 (1985). 

[4] K. Leung and D. Chandler, J. Stat. Phys. 63, 837 (1991). 

[5] A. P. Chatterjee, J. Chem. Phys. 113, 9310 (2000). 

[6] R. H. J. Otten and P. van der School, Phys. Rev. Lett. 
103, 225704 (2009). 

[7] R. H. J. Otten and P. van der School, J. Chem. Phys. 
134, 094902 (2011). 

[8] Z. Neda, R. Florian, and Y. Brechet, Phys. Rev. E 59, 
3717 (1999). 

[9] M. Foygel, R. D. Morris, D. Anez, S. French, and V. L. 
Sobolev, Phys. Rev. B 71, 104201 (2005). 

[10] L. Berhan and A. M. Sastry, Phys. Rev. E 75, 041120 
(2007). 

[11] T. Schilling, S. Jungblut, and M. A. Miller, Phys. Rev. 
Lett. 98, 108303 (2007). 

[12] G. Ambrosetti, C. Grimaldi, I. Balberg, T. Maeder, A. 
Danani, and P. Ryser, Phys. Rev. B 81, 155434 (2010). 

[13] R. M. Mutiso, M. G. Sherrott, J. Li, and K. 1. Winey, 
Phys. Rev. B 86, 214306 (2012). 

[14] B. Nigro, G. Grimaldi, P. Ryser, A. P. Ghatterjee, and 
P. van der School, Phys. Rev. Lett. 110, 015701 (2013). 

[15] A. P. Ghatterjee, J. Chem. Phys. 132, 224905 (2010); J. 
Stat. Phys. 146, 244, (2012). 

[16] 1. Balberg, C. H. Anderson, S. Alexander, and N. Wag¬ 
ner, Phys. Rev. B 30, 3933 (1984). 

[17] T. Schilling, M. Miller, and P. van der School, 
arXiv: 1505.07660; H. Meyer, P. van der School, and T. 


Schilling, J. Chem. Phys. 143, 044901 (2015). 

[18] M. Penrose, Random Geometric Graphs (Oxford Univer¬ 
sity Press, New York, 2003). 

[19] J. Dali and M. Christensen, Phys. Rev. E 66, 016121 

( 2002 ). 

[20] C. Herrmann, M. Barthelemy, and P. Provero, Phys. Rev. 
E 68, 026128 (2003). 

[21] N. Masuda, H. Miwa, and N. Konno, Phys. Rev. E 71, 
036108 (2005). 

[22] C. Grimaldi, Phys. Rev. E 92, 012126 (2015). 

[23] J. P. Straley, Phys. Rev. A 8, 2181 (1973). 

[24] L. Onsager, Ann. N.Y. Acad. Sci., 51, 627 (1949). 

[25] R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 

( 2002 ). 

[26] M. E. J. Newman, S. H. Strogatz, and D. J. Watts, Phys. 
Rev. E 64, 026118 (2001). 

[27] E. A. Leicht and R. M. DSouza, arXiv:0907.0894 

[28] R. Gohen, K. Erez, D. ben-Avraham, and S. Havlin Phys. 
Rev. Lett. 85, 4626 (2000). 

[29] B. Nigro, G. Grimaldi, M. A. Miller, P. Ryser, and T. 
Schillin g , Phys. Rev. E 88, 042140 (2013). 

[30] R. M. Mutiso and K.I. Winey, Prog. Polym. Sci. 40, 63 
(2015). 

[31] B. Vigolo, G. Goulon, M. Mangey, G. Zakri, and P. 
Poulin, Science 309, 920 (2005). 

[32] J.-P. Hansen and 1. R. McDonald, Theory of Simple Liq¬ 
uids (Elsevier, New York, 2006). 

[33] M. Franco-Melgar, A. J. Haslam, and G. Jackson, Mol. 
Phys. 106, 649 (2008). 

[34] A. P. Chatterjee, J. Chem. Phys. 141, 034903 (2014). 

[35] J. G. Meier, C. Crespo, J. L. Pelegay, P. Castell, R. 
Sainz, W. K. Maser, and A. M. Benito, Polymer 52, 1788 
( 2011 ). 



